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The Fat-Link Irrelevant Clover (FLIC) fermion action provides a new form of nonperturbative O(a) improve- 
ment and allows efficient access to the light quark-mass regime. FLIC fermions enable the construction of the 
nonperturbatively 0(a)-improved conserved vector current without the difficulties associated with the fine tuning 
of the improvement coefficients. The simulations are performed with an 0(a 2 ) mean-field improved plaquette- 
plus-rectangle gluon action on a 20 3 x 40 lattice with a lattice spacing of 0.128 fm, enabling the first simulation 
of baryon form factors at light quark masses on a large volume lattice. Magnetic moments, electric charge radii 
and magnetic radii are extracted from these form factors, and show interesting chiral nonanalytic behavior in the 
light quark mass regime. 



1. INTRODUCTION 

The magnetic moments of baryons have been 
identified as providing an excellent opportunity 
for the direct observation of chiral nonanalytic 
behavior in lattice QCD, even in the quenched 
approximation |1I2I3| . This paper will present 
results for baryon electromagnetic structure in 
which the chiral nonanalytic behaviour predicted 
by quenched chiral perturbation theory is ob- 
served in the numerical simulation results. 

The numerical simulations of the electromag- 
netic form factors presented here are carried 
out using the Fat Link Irrelevant Clover (FLIC) 
fermion action |4I5| in which the irrelevant oper- 
ators introduced to remove fermion doublers and 
lattice spacing artifacts JH| are constructed with 
smoothed links. These links are created via APE 
smearing [7J . On the other hand, the relevant op- 
erators surviving in the continuum limit are con- 
structed with the original untouched links gener- 
ated via standard Monte Carlo techniques. 

FLIC fermions provide a new form of non- 
perturbative Ola) improvement |5I8| where near- 
continuum results are obtained at finite lattice 
spacing. Access to the light quark mass regime 
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is enabled by the improved chiral properties of 
the lattice fermion action The 0(a)-improved 
conserved vector current [9f is used. Nonper- 
turbative improvement is achieved via the FLIC 
procedure where the terms of the Noether cur- 
rent having their origin in the irrelevant opera- 
tors of the fermion action are constructed with 
mean-field improved APE smeared links. The 
use of links in which short-distance fluctuations 
have been removed simplifies the determination of 
the coefficients of the improvement terms in both 
the action and its associated conserved vector 
current. Perturbative renormalizations are small 
for smeared links and are accurately accounted 
for by small mean-field improvement corrections. 
Hence, we are able to determine the form factors 
of octet and decuplet baryons with unprecedented 
accuracy. 

2. LATTICE TECHNIQUES 

2.1. Gauge and Quark Actions 

The simulations are performed using a mean- 
field C(a 2 )-improved Luscher-Weisz ^U] gauge 
action on a 20 3 x 40 lattice with a lattice spacing 
of 0.128 fm as determined by the Sommer scale 
ro = 0.50 fm. We use a minimum of 255 config- 
urations and the error analysis is performed by a 
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third-order, single-elimination jackknife. 

For the quark fields, we use the Fat-Link Ir- 
relevant Clover fermion action 0j. Fat links are 
created using APE smearing \7\. The smearing 
procedure replaces a link, U^(x), with a sum of 
the link and a times its staples 
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followed by projection back to SU(3). We select 
the unitary matrix J7j L which maximizes 

by iterating over the three diagonal SU(2) sub- 
groups of SU(3). We repeat this procedure of 
smearing followed immediately by projection n 
times. We select a smearing fraction of a — 0.7 
(keeping 0.3 of the original link) and iterate the 
process six times The mean- field improved 

FLIC 4- action now becomes 

Sik = S& ~ tWa^F^x) , (2) 

where F^ v is constructed using fat links, and 
where the mean-field improved Fat-Link Irrele- 
vant Wilson action is 



(3) 



7/x 



"0 



ij)(x + A) 



utix - A) 



ip(x ~ A) 



ip(x + A) 



U^(x-H) 



,FL 



^(a; - A) 



with k = l/(2m + 8r). We take the standard 
value r = 1. The 7-matrices are hermitian and 
0- M „ = [7^, 7^]/(2«)- 



For fat links, the mean link Uq rj 1, enabling 
the use of highly improved definitions of the lat- 
tice field strength tensor, ,6J- In particular, 
we employ an 0(a 4 )-improved definition of F^ 
in which the standard clover-sum of four lxl 
Wilson loops lying in the fi, v plane is combined 
with 2x2 and 3x3 Wilson loop clovers. More- 
over, mean-field improvement of the coefficients 
of the clover and Wilson terms of the fermion ac- 
tion is sufficient to accurately match these terms 
and eliminate O(a) errors from the fermion action 

.Previous work PiTHl has shown that the FLIC 
fermion action has extremely impressive conver- 
gence rates for matrix inversion, which provides 
great promise for performing cost effective sim- 
ulations at quark masses closer to the physical 
values. Problems with exceptional configurations 
have prevented such simulations in the past. The 
ease with which one can invert the fermion matrix 
using FLIC fermions leads us to attempt simula- 
tions down to light quark masses corresponding 
to m^/nip = 0.35 in an attempt to reveal chiral 
non-analytic behaviour in baryon magnetic mo- 
ments. 

A fixed boundary condition at t = is used 
for the fermions and gauge-invariant Gaussian 
smearing |13I14| in the spatial dimensions is ap- 
plied at the source at t = 8 to increase the over- 
lap of the interpolating operators with the ground 
state while suppressing excited state contribu- 
tions. 

2.2. Improved Conserved Vector Current 

For the construction of the C(a)-improved con- 
served vector current, we follow the technique 
proposed by Martinelli et al. 23 . The standard 
conserved vector current for Wilson-type fermions 
is derived via the Noether procedure 



3 



p. r)Up(x)rp(x + A) 

+ ^(x + fx)(^ + r)Ul(x)i)(x) 

+ (x^x-fi)}. (4) 

The 0(a) improvement term is also derived from 
the fermion action and is constructed in the form 
of a total four-divergence, preserving charge con- 
servation. The 0(a)-improved conserved vector 
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current is 
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j^( x ) + ^ c cvc a^d p (ilj(x)a p ^{x)) ,(5) 



where Ccvc is the improvement coefScient for 
the conserved vector current and we define 

d p @(x)i/>(x)) = ?(x)(^ p + ^ P )^(x) , (6) 

where the forward and backward derivatives are 
defined as 
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The terms proportional to the Wilson parame- 
ter r in Eq. Q and the four-divergence in Eq. J5J 
have their origin in the irrelevant operators of 
the fermion action and vanish in the continuum 
limit. Nonperturbative improvement is achieved 
by constructing these terms with fat-links. As 
we have stated, perturbative corrections are small 
for fat-links and the use of the tree-level value 
for Ccvc — 1 together with small mean-field im- 
provement corrections ensures that O(a) artifacts 
are accurately removed from the vector current. 
This is only possible when the current is con- 
structed with fat-links. Otherwise, Ccvc needs 
to be appropriately tuned to ensure all 0(a) ar- 
tifacts are removed. 

2.3. Lattice Three-Point Functions 

The technique used for constructing the three- 
point functions follows the procedure outlined in 
detail in Refs. |15I16| . In particular, we use the 
sequential source technique at the current inser- 
tion. Correlation functions are made purely real 
with exact parity through the consideration of U 
and U* link configurations. Electric and magnetic 
form factors are extracted by constructing ratios 
of two- and three-point functions for a baryon, B 

\c B ^ B (t 2M -,P ,m){c B ^ B (t 2 M-,-P,-P'-x))\ 1/2 

(G BB (t 2 ;P;T i ))(G BB (t 2 ;-p-r i )) / 

-(^) 1/2 (^) 1/2 ^,p;r;M) (7) 



where > ti > <2 are the current insertion and 
baryon sink time slices respectively; p and p' are 
the baryon momentum before and after current 
insertion respectively; and T is a 4 x 4 dirac ma- 
trix with T4 = |(l+74)- The Sachs form of the 
electromagnetic form factors, 
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are obtained from R by 
R(q,6;UA) = Q E {q 2 ) 



Guiiq 2 ) Vijkq 1 
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R(q,0;T v k) = 

*®ff;r<,*) = {Eq + M) . 

We simulate at the smallest finite q 2 available on 
our lattice, q = ^j-x. We insert our improved con- 
served vector current at t\ — 14. Since R is inde- 
pendent of the baryon sink time slice, we calculate 
the ratio in Eq. Q for a range of sink times. Us- 
ing a covariance- matrix-based \ 2 /Ndof, we inde- 
pendently select suitable fitting windows for both 
electric and magnetic form factors. 

2.4. Correlation Function Analysis 

In selecting the most appropriate Euclidean 
time window for fitting, we consider the corre- 
lated x 2 /dof as given by 
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where, M (= 1) is the number of parameters to be 
fitted, Nt is the number of time slices considered, 
y(ti) is the configuration average value of the de- 
pendent variable at time tj that is being fitted to 
a constant value k, and C(ti, tj) is the covariance 
matrix. The elements of the covariance matrix 
are estimated via the jackknife method, 
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Table 1 

Fit results for the contribution of a u quark of 
unit charge to the proton magnetic form factor 
(/ijy) obtained for various time-fitting windows. 



Time 


Fit Value 


Uncertainty 


X 2 /dof 


16-18 


0.865 


0.022 


7.891 


16-19 


0.867 


0.024 


6.304 


16-20 


0.868 


0.025 


4.731 


17-19 


0.884 


0.03 


1.281 


17-20 


0.885 


0.033 


0.855 


18-20 


0.893 


0.042 


0.047 




where, N c is the total number of configurations, 
Vmifi) is the jackknife ensemble average of the 
system after removing the m t h configuration. 
y(U) is the average of all such jackknife averages, 
given by 

1 Nc 

m— 1 

Tabled shows the contribution from a u quark 
of unit charge to the proton magnetic form fac- 
tor in nuclear magnetons (/it at). The uncertainty 
and the x 2 / dof are quoted for various time-fitting 
windows. The principal value of the fitted param- 
eter remains fairly constant over different time 
windows, but the x 2 /dof shows a marked de- 
crease as we move from the time window 16-18 
to 17-20. Figure ^ illustrates these lattice simula- 
tion results obtained from Eq. (JJJ) for one of the 
lighter quark masses. The horizontal line indi- 
cates the best fit in the time window 17-20 which 
provides the most acceptable x 2 /dof ~ 1. 

3. RESULTS 

Figure displays FLIC fcrmion simulation re- 
sults for the magnetic moments of the proton and 
neutron in quenched QCD. The moments are ob- 
tained using the empirical fact that the Q 2 de- 
pendence of the electric and magnetic form fac- 
tors of the proton over the range — > 0.23 GeV 2 
are nearly equivalent. At heavy quark masses we 
note that the magnetic moments of both nucle- 
ons display approximately linear behaviour when 



Figure 1. Contribution from a u quark of unit 
charge to the proton magnetic form factor. The 
quark mass for this correlation functions corre- 
sponds to ml ~ 0.3 GeV 2 . 



plotted as a function of m\. Simulations at mod- 
erately heavier quark masses are expected to re- 
veal a Dirac moment behavior oc l/m q ~ 1/m 2 ,. 
As one approaches the light quark mass regime, 
we find evidence of non-analytic behaviour in 
the nucleon magnetic moments as predicted by 
quenched %PT |1I3| . In fact, if we were to flip the 
sign of the neutron magnetic moment, we would 
find that the proton and neutron have a very sim- 
ilar behaviour as a function of m\ in the light 
quark mass regime as predicted by the leading 
non-analytic contributions of quenched xPT. 

Figures E3 and 0] display results for the proton 
and neutron charge radii, respectively, obtained 
from a dipole form factor ansatz. For the proton 
some curvature is emerging as the chiral limit is 
approached. The FLIC fermion simulations re- 
veal a small, negative value for the calculation 
of the quenched neutron charge radius at large 
and intermediate quark masses before the signal 
is lost at the lightest two quark masses. This re- 
sult confirms the earlier result that the two 
d quarks have a larger charge radius than the u 
quark within the neutron. This is also in agree- 
ment with quark model calculations |17j . 

We also perform a calculation of the individual 
quark sector contributions to the magnetic mo- 
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Figure 2. FLIC fermion simulation results for the 
magnetic moments of the proton (upper) and neu- 
tron (lower) in quenched QCD. The experimental 
values are given by asterisks 
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Figure 3. FLIC fermion simulation results for the 
charge radius of the proton in quenched QCD. 



mcnts. In Fig. [5] we show results for a calculation 
in quenched QCD of the ratio of singly to doubly 
represented quark magnetic contributions in the 
nucleon for quarks of unit charge. We see that 
our results deviate significantly from the simple 
quark model prediction of —1/2. 

Figure [5] shows the results for a calculation of 
the octet S baryons. We note the level order- 
ing of the three E baryons, with the positive and 
negative charge states both approaching the ex- 
perimental values at light quark masses. We also 
see that our quenched lattice calculation of the 
magnetic moment of the neutral E baryon pre- 
dicts a value in the range 0.5 — 0.7fiN, although 
we first need to take into account the appro- 
priate chiral extrapolation and quenching effects 
|18j before an accurate prediction can be made. 
Our results reveal a significant amount of nonana- 
lytic behaviour for the E + , which is predicted by 
quenched xPT 0- QxPT also predicts a small 
amount of curvature for the E° and E~ and we 
confirm this prediction 3 . 
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Figure 4. FLIC fermion simulation results for the 
charge radius of the neutron in quenched QCD. 
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Figure 5. Ratio of singularly to doubly repre- 
sented quark magnetic contributions in the nu- 
cleon for quarks of unit charge. 
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Figure 7. FLIC fermion simulation results for 
the magnetic moment of the octet 5 baryons in 
quenched QCD. 
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Figure 6. FLIC fermion simulation results for 
the magnetic moment of the octet £ baryons in 
quenched QCD. 



We have also performed a calculation of the 
magnetic moment of the octet S baryons in 
quenched QCD and these results are shown in 
Fig. \7\ Here we see an improved signal for the 
magnetic moment of the 3 compared to the other 
baryons at light quark masses due to the presence 
of two strange quarks. We note the splitting be- 
tween the neutral and negative charged baryons is 
close to the experimentally observed value, even 
in the quenched approximation. Predictions from 
quenched chiral perturbation theory P] suggest 
that there should be no nonanalytic behaviour 
for the S~ and a small amount of downward cur- 
vature for the 2°. Our results confirm these pre- 
dictions. 

Figure |H1 displays FLIC fermion simulation re- 
sults for the magnetic moments of the proton and 
A + resonance in quenched QCD. At large pion 
masses, the A moment is enhanced relative to the 
proton moment in accord with earlier lattice QCD 
results |15llfi| and model expectations. However 
as the chiral regime is approached the nonanalytic 
behavior of the quenched meson cloud is revealed, 
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Figure 8. FLIC fermion simulation results for the 
magnetic moments of the proton (o) and A + res- 
onance (A) in quenched QCD. 



enhancing the proton and suppressing the A + in 
accord with the expectations of quenched xPT 
|19I20| . The quenched artifacts of the A provide 
an unmistakable signal for the onset of quenched 
chiral nonanalytic behavior. 

4. CONCLUSIONS 

We have presented the first lattice QCD simu- 
lation results for the electromagnetic form fac- 
tors of the nucleon, E, 3 and A baryons at 
quark masses light enough to reveal unmistakable 
quenched chiral nonanalytic behavior. The non- 
linear behaviour observed is in agreement with 
the predictions of quenched chiral perturbation 
theory |5I20] . 

This work paves the way for a study of the indi- 
vidual quark sector contributions to baryon mag- 
netic moments and, in particular, the strangeness 
content of the nucleon J5| . 

This work was supported by the Australian Re- 
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generous grants of supercomputer time which 
have enabled this project. 

REFERENCES 

1. D. B. Leinweber, Nucl. Phys. Proc. Suppl. 
109 (2002) 45 [hep-lat/01 12021) . 

2. M. J. Savage, Nu cl. Phys. A 700 (2002) 359 
nucl-th/0107038 . 

3. D. B. Leinweber, to appear in Phys. Rev. D, 
hep-lat/0211017. 

4. J.M. Zanotti et al, Phys. Rev. D 
60 (2002) 074507 |hep-lat/0110216| ; 
Nucl.Phys.Proc.S uppl. 109 101 (2002) 
|hep-lat/020i004) . 

5. D. B. Leinweber, et al, |nucl-th/0211014| 

6. S. O. Bilson-Thompson, D. B. Leinweber and 
A. G. Williams, Annals Phys. 304 (2003) 1 
[hep-lat/0 203008 . 

7. M. Falcioni et al, Nucl. Phys. B251 (1985) 
624; M. Albanese et al, Phys. Lett. B 192 
(1987) 163. 

8. J. Zanotti, et al, in preparation. 

9. G. Martinelli, C. T. Sachrajda and 
A. Vladikas, Nucl. Phys. B 358 (1991) 
212. 

10. M. Luscher and P. Weisz, Commun. Math. 
Phys. 97, 59 (1985) [ibid. 98, 433 (1985)]. 

11. F.D. B onnet et al, Phys . Rev. D 62 (2000) 
094509 hep-lat/0001018 . 

12. W. Kamleh, D. H. Adams, D. B. Leinweber 
and A. G. Williams, Phys. Rev. D 66, 014501 

(2002) |hep-lat/011204l) . 

13. S. Gusken, Nucl. Phys. Proc. Suppl. 17, 361 
(1990). 

14. J. M. Zanotti, et al, Phys. Rev. D 68, 054506 

(2003) |hep-lat/030400ll . 

15. D. B. Leinweber, R. M. Woloshyn and 
T. Draper, Phys. Rev. D 43 (1991) 1659. 

16. D. B. Lein weber et al, Phy s. Rev. D 46 
(1992) 3067 hep-lat/9208025 . 

17. R. D. Carlitz, S. D. Ellis and R. Savit, Phys. 
Lett. B 68, 443 (1977). 

18. D. B. Leinweber, et al, these proceedings. 

19. D. B. Leinweber, et al, nucl-th/0308083. 

20. R. D. Young, D. B. Leinweber and 
A. W. Thoma s, |hep-lat/0309187| and 
|hep-lat/0311038 



